An Improved Exact Algorithm for the 
Domatic Number Problem* 



Tobias Riege, ^ Jorg Rothe, ^ Holger Spakowski ^ 
Institut fur Informatik 
Heinrich- Heine- Universitdt Diisseldorf 
Diisseldorf, Germany 

Masaki Yamamoto ^ 
Department of Mathematical and Computing Science 
Tokyo Institute of Technology 
Tokyo, Japan 

March 16, 2006 



Abstract 

The S-domatic number problem asks whether a given graph can be partitioned into 
three dominating sets. We prove that this problem can be solved by a deterministic 
algorithm in time 2.695" (up to polynomial factors). This result improves the previous 
bound of 2.8805"', which is due to Fomin, Grandoni, Pyatkin, and Stepanov. To prove 
our result, we combine an algorithm by Fomin et al. with Yamamoto's algorithm for the 
satisfiability problem. In addition, we show that the 3-domatic number problem can be 
solved for graphs G with bounded maximum degree A(G') by a randomized algorithm, 
whose running time is better than the previous bound due to Riege and Rothe [RR05] 
whenever A(G) > 5. Our new randomized algorithm employs Schoning's approach to 
constraint satisfaction problems. 
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1 Introduction 



A dominating set in an (undirected) graph G is a subset D of the vertex set F of G such 
that the closed neighborhood of D equals V. The domatic number problem asks for a 
partition of V into a maximum number of disjoint dominating sets. This number, denoted 
by S{G), is called the domatic number of G. The domatic number problem arises in the area 
of computer networks and is related to the tasks of distributing resources and of locating 
facilities in the network. This problem has been intensely studied, see, e.g., [CH77, Far84, 
Bon85, KS94, HT98, FHKOO, FHKS02, RR04]. 

For each /c > 3, it is NP-complete to determine whether or not the domatic number of 
a given graph is at least fe, see Garey and Johnson [GJ79]. That is why we cannot expect 
to find a polynomial-time algorithm that solves the problem. 

Among the various ways of coping with NP-hard problems (such as approximation, ran- 
domization, or parameterized complexity), much attention has been paid to designing exact 
exponential-time algorithms for such problems that are better than the trivial exponential- 
time algorithm. In particular, if the trivial algorithm runs in time 3" but one is able to 
find a algorithm for this problem with c < 3 (up to polynomial factors), then this al- 
gorithm can handle larger problem instances in the same amount of time than the trivial 
algorithm. This difference can be quite important in practice. For recent surveys on this 
subject, we refer to Schoning [Sch05] and Woeginger [Woe03]. 

The first result breaking the trivial 3" barrier for the 3-domatic number problem is due to 
Ricgc and Rothe [RR05]. Fomin et al. [FGPS05] improved their result by providing a 2.8805" 
bound. In this paper, we obtain a further improvement by making use of two known 
algorithms: an algorithm by Fomin et al. [FGPS05] for generating all minimal dominating 
sets of a graph, and Yamamoto's algorithm for the satisfiability problem [Yam05]. We 
show that these two algorithms can be combined so as to yield a 2.695" time algorithm for 
the 3-domatic number problem. 

In addition, we prove that there is a randomized algorithm solving this problem for 
graphs with bounded maximum degree, which improves the previous bound due to Riege 
and Rothe [RR05]. Here, we apply Schoning's results for constraint satisfaction prob- 
lems [Sch99, Sch02], which previously was also useful in improving the bounds of random- 
ized exponential-time algorithms for the satisfiability problem. For example, Iwama and 
Tamaki [IT04] designed a randomized algorithm with running time 0(1.324") for solving 
3-SAT by combining Schoning's algorithm with an algorithm due to Paturi et al. [PPSZ98]. 

2 Preliminaries 

We first introduce some graph-theoretical notation. Graphs are pairs G = {V,E), where V 
is the vertex set of G and E is the edge set of G. All graphs considered in this paper are 
undirected and simple. That is, the edges of any graph are unordered pairs of vertices, 
and there are neither multiple nor reflexive edges; i.e., there exists at most one edge of the 
form {x, y} for any two vertices x and y, and there is no edge of the form {x, x^. Moreover, 
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we require all graphs to not have isolated vertices. In general, graphs need not be connected. 
Definition 1 Let G = {V, E) be a graph. 

• For any vertex v eV, define the neighborhood of f in G by 

N{v) = {ueV\{u,v} e E}, 
and define the closed neighborhood of u in G by 

N[v] = {v}UN{v). 

• For any subset U C. V of the vertices of G, define N\U] = U«g!7-^M '^''^'^ 
N{U) = N[U] - U. 

• A subset D C. V is a dominating set of G if and only if every vertex u ^ D is a 
neighbor of some vertex v E: D. That is, D dominates G if and only if N[D] = V. 

• A dominating set D is called a minimal dominating set if and only if there exists no 
dominating set C of G with C C D. 

• The domatic number of G (denoted by S{G)) is the maximum number of disjoint 
dominating sets. 

• Given a graph G and a positive integer k, the fc-domatic number problem, k-DNP for 
short, asks whether or not S{G) > k. 

Note that at least one partition of G into S{G) (i.e. the maximum number of) disjoint 
dominating sets contains 6{G) — 1 minimal dominating sets. 

For A; > 3, the A;-domatic number problem is NP-complete, see Garey and John- 
son [GJ79]. Therefore, no polynomial-time algorithm for k-DNP exists unless P = NP. 

Here, we focus on the case k = 3 only. The first algorithm beating the trivial 0(3") 
barrier^ for the 3-domatic number problem is due to Riege and Rothe [RR05]. They also 
investigate this problem for graphs with bounded maximum degree, and they propose a 
deterministic and a randomized algorithm for it. In particular, this deterministic algo- 
rithm outperforms their general deterministic algorithm for the 3-domatic number problem 
whenever A(G), the maximum degree of the input graph G, is in the range 3 < A(G) < 7. 

Fomin et al. [FGPS05] recently improved the general result from [RR05] by constructing 
an algorithm for 3-DNP with running time 0(2.8805""). Their algorithm makes use of a new 
algorithm for enumerating all minimal dominating sets of a graph. 

Theorem 2 (Fomin et al. [FGPS05]) There is an algorithm for listing all minimal 
dominating sets in an n vertex graph G in time 0(1.7697"). 

^As is common for exponential-time algorithms, we use the O notation to indicate that polynomial factors 
are neglected. That is, for functions / and g, we write / £ if / € 0{p ■ g), where p is some polynomial. 
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Their proof relies on a new method to evaluate the size of the recursion tree for a given 
exponential-time algorithm. This technique, which is called "measure and conquer," was 
introduced by Fomin, Grandoni, and Kratsch [FGK05], who applied it to give a better 
analysis of the runtimes of exact backtracking algorithms for the minimum dominating set 
problem and the minimum set cover problem. 

The "measure and conquer" technique is based on the following idea. By choosing a 
suitable measure of the subproblems generated by the recursive algorithm considered, one 
can lower-bound the progress made by the algorithm in each branching step. A clever choice 
of this measure can yield a much better worst-case runtime analysis of the problem, even 
though the algorithm considered is not new and has long been known. 

Based on this new technique, Fomin et al. [FGPS05] design an cxponcntial-time al- 
gorithm for determining the domatic number of a given graph (and thus, in particular, 
for solving 3-DNP). Their approach resembles the dynamic-programming algorithm by 
Lawler [Law76] for computing the chromatic number of a graph. 

Corollary 3 (Fomin et al. [FGPS05]) There is an algorithm for computing the domatic 
number of a given graph that runs in time 0(2.8805"). 

Prior to this paper, the time bound stated in Corollary 3 was the best result known for 
the (search version of the) domatic number problem, and in particular for 3-DNP. In the 
next section, we further improve the result that follows from Corollary 3 for the decision 
problem 3-DNP. 

3 An Improved Exact Algorithm for the 3-Domatic Number 
Problem 

In this section, we improve the exponential running time for 3-DNP that follows from Corol- 
lary 3. To this end, we combine the algorithm for enumerating all minimal dominating sets 
from Theorem 2 with an algorithm due to Yamamoto [Yam05]. 

Yamamoto's algorithm, which is based on (and improves) an algorithm due to 
Hirsch [HirOO], solves the NP-complete problem SAT in time 0(1.234™), where m is the 
number of clauses of the given boolean formula in conjunctive normal form. To prove The- 
orem 5 below, we will apply Yamamoto's algorithm to a special version of SAT, which is 
called NAE-SAT ("not-all-equal satisfiability"). 

Definition 4 Let if = ip{X, C) he a boolean formula in conjunctive normal form consisting 
of a collection C = {ci, C2, . . . , c^} of m clauses over the variable set X. We say is 
in NAE-SAT if and only if there exists a truth assignment for X satisfying all clauses in C 
and such that in none of the clauses, all literals are true. 

Theorem 5 There exists an exact algorithm solving the 3-DNP problem in time 0(2.695"). 
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Proof. Wc will only sketch the proof. 

Let G = iy, E) be a given graph with n vertices. Using the algorithm from Theorem 2, 
generate all minimal dominating sets D C. V ol G. Given one such minimal dominating 
set -D, create a formula ipo = 'Pd{X,G) for the NP-compIete problem NAE-SAT as follows: 

• The set of variables is defined as 

X = {xy\ve V}. 

• For each vertex v eV, create the clause 



[J Xu> , 

u e N[v\ 

so the clause set is defined as C = {G^ \v ^V}. 

Note that G G 3-DNP if and only if (^d is in NAE-SAT for some formula tpo thus defined. 

(Alternatively, one might formalize this 3-DNP instance G with some minimal dominating 
set D as an instance {U,Sd) of the minimum set cover problem, where U = {0, 1,2} is a 
universe of elements (which correspond to the three dominating sets of G if G is in 3-DNP) , 
and where Sd is a correspondingly defined nonempty set system over U. Then, we assign 
the element 2 of U to D, and G G 3-DNP if and only if every set S in Sd has both the 
values and 1. This corresponds exactly to the NAE-SAT property stated above.) 

The number of clauses in equals n, the number of vertices in G. Using the standard 
polynomial-time reduction from NAE-SAT to SAT (see Garey and Johnson [GJ79]), which 
adds the negation of each clause to the clause set, we obtain a formula (p'j-, = {X, C U C) 
with the property that 

(PD & NAE-SAT Lf'j^ G SAT. 

Using the exponential-time SAT algorithm designed by Yamamoto [Yam05], we can now 
determine the satisfiability of ip'j^ in time 0(1. 234^"). Note here that tp'jy has 2n clauses. It 
follows that our algorithm to solve 3-DNP has a running time of 

0(1.7697" • 1.234^") = 0(2.695"), 

which completes the proof. | 



4 An Improved Randomized Algorithm for the 3-Domatic 
Number Problem 

We now turn to the case where the maximum degree A(G) of the input graph G is bounded 
by some small constant. Riege and Rothe [RR05] described a randomized algorithm for this 
problem. Here, we observe that their result can be improved by using Schoning's algorithm 
for constraint satisfaction problems [Sch99, Sch02]. 
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A(G) 


3 


4 


5 


6 


7 


8 


source 


deterministic 


2.2894" 


2.6591" 


2.8252" 


2.9058" 


2.9473" 


2.9697" 


[RR05] 


randomized 


2" 


2.3570" 


2.5820" 


2.7262" 


2.8197" 


2.8808" 


[RR05] 


randomized 


2.2501" 


2.4001" 


2.5001" 


2.5715" 


2.6251" 


2.6667" 


Thm. 6 



Table 1: Three 3-DNP algorithms for graphs G with bounded maximum degree A(G). 



Theorem 6 There is a randomized algorithm solving the 3-DNP problem for graphs G with 
bounded maximum degree A(G), whose running time depends on A{G) as stated in Table 1. 

Proof. The 3-DNP problem for a graph G can easily be formulated as a constraint sat- 
isfaction problem (CSP) on the domain D = {0,1,2} with each constraint having or- 
der £ = A(G) -|- 1. The order of a constraint is the number of arguments in the constraint. 

Let G = (y, E) be a given graph. Create the constraint satisfaction problem F with n 
variables as follows 

• The set of variables is defined as 

X = {xy\ve V}. 

• For each vertex -y G F, create the constraint defined by 

Gv[Xy, Xu)l , Xw2 1 ■ ■ ■ ■> ^«'||JV[!,]||-l) ~ -'■ 

all values oi D = {0,1,2} appear in the values of Xy^Xy]^^Xyj2i ■ ■ ■ ^ Xiu,,^,-,,,^., 

where iwi, ■«;2, • • • , f«||iv[i;]||-i the vertices adjacent to v. 
It is easy to see that 

G G 3-DNP <J=^ CSP F has a solution. 

To determine whether F has a solution, we can apply Schoning's randomized algorithm 
for solving constraint satisfaction problems [Sch99, Sch02]. 

We thus obtain a randomized algorithm for 3-DNP with running time 



O 3 1 



A(G) + 1 



+ e 



for any e > 0. I 

For graphs with A(G) > 5, this is an improvement over the randomized algorithm given 
by Riege and Rothe [RR05], see Table 1. 
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5 Conclusions 



In this paper, we considered the 3-domatic number problem, which asks whether a given 
graph can be partitioned into three dominating sets, and we have shown how to improve 
on existing exact and randomized exponential-time algorithms for this problem. 

In particular, by reducing 3-DNP to the NAE-SAT problem and by combining Yamamoto's 
algorithm [Yam05] with the algorithm by Fomin et al. [FGPS05], we obtained an exact (i.e., 
deterministic) algorithm that runs in time 0(2.695"). This result improves on the previously 
best bound of 0(2.8805"') for the 3-domatic number problem, which was achieved by Fomin 
et al. [FGPS05]. An even earlier — and to our knowledge, the first — nontrivial bound for 
this problem is due to Riege and Rothe [RR05] who presented an 0(2.9416'^) algorithm. 

We also described a randomized algorithm solving the 3-domatic number problem for 
graphs G with bounded maximum degree A(G). Whenever A(G) > 5, the running time of 
our new randomized algorithm is better than the previously known bound, which is due to 
Riege and Rothe [RR05]. These results are summarized in Table 1. Our new randomized 
algorithm makes use of Schoning's algorithm for constraint satisfaction problems [Sch99, 
Sch02]. 
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ing discussions on the subject of this paper. In particular, we thank Dieter Kratsch for 
calling our attention to Lawler's algorithm. The second author thanks Osamu Watanabe 
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